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Abstract 

Previous  studies  have  shown  that  the  amplitude  of  the  vibration  of  a  piezoelectric  (PZT)  device  produces  an  oscillating  flow  that  changes  the 
chamber  volume  along  with  a  curvature  variation  of  the  diaphragm.  In  this  study,  an  actuating  micro-diaphragm  with  piezoelectric  effects  is  utilized 
as  an  air-flow  channel  in  proton  exchange  membrane  fuel  cell  (PEMFC)  systems,  called  PZT-PEMFC.  This  newly  designed  gas  pump,  with  a 
piezoelectric  actuation  structure,  can  feed  air  into  the  system  of  an  air-breathing  PEMFC.  When  the  actuator  moves  outward  to  increase  the  cathode 
channel  volume,  the  air  is  sucked  into  the  chamber;  moving  inward  decreases  the  channel’s  volume  and  thereby  compresses  air  into  the  catalyst  layer 
and  enhancing  the  chemical  reaction.  The  air- standard  PZT-PEMFC  cycle  is  proposed  to  describe  an  air-breathing  PZT-PEMFC.  A  novel  design 
for  PZT-PEMFCs  has  been  proposed  and  a  three-dimensional,  transitional  model  has  been  successfully  built  to  account  for  its  major  phenomena 
and  performance.  Moreover,  at  high  frequencies,  PZT  actuation  leads  to  a  more  stable  current  output,  more  drained  water,  higher  sucked  air,  higher 
hydrogen  consumption,  and  also  overcomes  concentration  losses. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Previous  studies  have  shown  that  air-breathing  proton 
exchange  membrane  fuel  cells  (PEMFC s)  are  lighter  in  weight 
because  they  do  not  require  an  additional  oxygen  feeding  device. 
In  addition,  the  performance  of  air-breathing  PEMFCs,  due  to 
oxygen  transport  limitations,  is  proportional  to  the  oxygen  trans¬ 
fer  coefficient  [1].  Santa  Rosa  et  al.  [2]  developed  an  open-air 
cathode  PEMFC  stack  with  forced  air  convection  instead  of 
natural  convection.  The  stack  performance  had  an  optimum  per¬ 
formance  for  the  air  fan  voltage  supply  at  5.0  V  when  using  an 
air  fan  for  combined  oxidant  supply  and  stack  cooling.  Zhang  et 
al.  [3]  presented  a  comprehensive  numerical  analysis  on  the  var¬ 
ious  stack  and  cell  parameters  for  stack  designs  with  an  array  of 
air-breathing  fuel  cell  cartridges.  Matamoros  and  Bruggemann 
[4]  showed  that  dehydrating  phenomena  slightly  affected  the 
performance  of  air-breathing  PEMFCs. 
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In  general,  catalyst  loading,  relative  humidity,  temperature, 
hydrogen  stoichiometry,  gas  diffusion  layer  (GDL)  thickness, 
and  cathode  structure  are  all  important  factors  in  the  perfor¬ 
mance  of  air-breathing  PEMFCs  [5,6].  Previous  studies  [7,8] 
have  indicated  that  interdigitated  flow  fields  force  the  reactant 
gas  to  pass  the  catalyst  layer,  which  increases  the  rate  of  chemi¬ 
cal  reaction  and  the  performance  of  fuel  cells,  thereby  reducing 
the  flooding  phenomena.  Thus,  in  this  study  the  interdigitated 
flow  field  is  applied  as  the  anode  channel.  Li  et  al.  [9]  pre¬ 
sented  the  design  of  a  flow  field  in  a  bi-polar  plate,  which  is  one 
of  the  key  components  in  PEMFC  stacks;  it  performs  a  num¬ 
ber  of  essential  functions  in  stack  operation,  such  as  supplying 
reactants  to  the  cell-active  area,  current  collection,  mechanical 
support  to  the  membrane  electrode  assembly  (MEA),  and  water 
and  heat  management,  as  well  as  the  maintenance  of  the  separate 
reactants. 

Micro-pumps  have  been  developed  with  several  actuation 
methods,  such  as  piezoelectric,  electromagnetic,  shape  mem¬ 
ory  alloy,  electrostatic,  and  thermo-pneumatic  devices  [10-18]. 
Most  have  complex  structures  and  low  pumping  rates.  However, 
piezoelectric  actuation  has  the  advantages  of  having  a  relatively 
simple  structure  and  low  power  consumption.  Micro-diaphragm 
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Nomenclature 

a't  stoichiometric  coefficient  of  the  products 

a!(  stoichiometric  coefficient  of  the  reactants 

Ain  inlet  area  (m2) 

Apzt  Piezoelectric  area  (m2) 

[A/y]eff  effective  surface  to  volume  ratio  (m-1) 

Deff  effective  diffusion  coefficient  (ms-1) 

/  frequency  of  PZT  (Hz) 

F  Faraday  constant  (96,439  C  mole-1) 
h  enthalpy  (KJ) 

i  current  (A) 

jo  reference  current  density  at  a  known  open  circuit 
voltage  (Am-3) 
jt  transfer  current  (A) 

Ji  species  diffusive  flux  (mol  m-2  s) 

ke ff  effective  thermal  conductivity  (W  m-1  K-1) 

M  Mass  (kg) 

N  amount  of  species 

P0  atmosphere  pressure  (N  m-2) 

R  gas  constant  ( J  mol- 1  K- 1 ) 

t  time  (s) 

T  temperature  (K) 

Va  anode  inlet  velocity  (ms-1) 

Vc  cathode  inlet  velocity  (ms-1) 

Vpzt  motion  equation  of  the  piezoelectric  device 
(ms-1) 

V  volume  (m3) 

Yi  mass  fraction  of  ith  species 

Greek  symbols 

a&,  ac  kinetic  constants  which  is  drawn  by  Tafel  curve 
through  the  experiments 
otk  concentration  exponents 

8  diffusion  length  (m) 

s  porosity 

rj  over  potential  (V) 

k  permeability  (m2) 

pc  channel  air  density  (kg  m-3) 

po  atmosphere  air  density  (kg  m-3) 

a  electrical  conductivity  (Q  m) 

r  shear  stress  (N  m-2) 

ip  potential  (V) 

(Pf  porous  (or  liquid)  potential  (V) 

<ps  solid  potential  (V) 

coi  production  rates  (kg  m-3  s-1 )  co  angular  velocity 
(rad  s-1) 

[A]  molar  concentration  of  reactant  in  porous  medium 

(mol) 


pumps  are  classified  based  on  the  presence  or  absence  of  valves 
[19,20].  A  plastic  micro-pump  with  a  valve  is  capable  of  pump¬ 
ing  both  liquid  and  gas  at  a  pumping  rate  of  2  ml  min-1  for 
water  and  up  to  50  ml  min-1  for  air  when  the  actuation  fre¬ 
quency  is  between  2  and  500  Hz  [21].  A  gas  pump  with  a 


novel  bimorph  actuation  structure  can  also  feed  air  into  the 
micro-direct  methanol  fuel  cell  (DMFC)  [22].  These  results 
showed  that  the  air-diaphragm  pump  worked  at  a  flow  rate  of 
85.3  ml  min-1  in  resonance  with  3.18  mW,  and  at  a  low-power 
consumption  of  below  20  V.  A  new  design  of  a  one-side  actuating 
micro-diaphragm  pump  with  PZT  was  successfully  developed 
that  exploited  the  harmonic  resonance  of  the  working  liquid 
with  the  other  system  components  (valve  and  diaphragm)  in 
the  pump  chamber  [23].  The  maximum  flow  rate  for  the  newly 
designed  pump  reached  72  ml  min-1  at  the  zero  pump  head  in 
the  operation  frequency  range  of  70-180  Hz  [24]. 

The  flow  channel,  gas  diffusion  layer,  catalyst  layer  and 
membrane  all  influence  PEMFC  performance.  During  PEMFC 
operation,  especially  at  low  stoichiometry,  liquid  water  is  likely 
to  appear  in  the  cathode,  reducing  oxygen  diffusion.  In  the 
membrane,  water  is  transported  by  electro-osmotic  drag,  back- 
diffusion,  and  convection.  Yi  and  Nguyen  [25]  showed  that  the 
performance  of  a  PEMFC  could  be  improved  by  anode  humid¬ 
ification  and  positive  differential  pressure  between  the  cathode 
and  anode,  to  increase  the  back  transport  rate  of  water  across 
the  membrane.  Ge  and  Yi  [26]  predicted  that  the  dry  reactant 
gases  can  be  successfully  humidified  internally  and  maintain 
high  performance  when  PEMFC  is  operated  in  the  counterflow 
mode  without  external  humidification.  The  electro-osmotic  drag 
coefficient  has  been  found  to  be  a  function  of  the  water  content  of 
the  membrane  and  temperature  [27-30].  Generally,  the  electro- 
osmotic  drag  coefficient  increases  with  increasing  water  content 
in  the  membrane.  At  high  water  content,  the  electro-osmotic 
drag  coefficient  increases  with  increasing  temperature.  Thus, 
water  and  thermal  management  is  a  key  to  successful  PEM- 
FCs,  especially  to  PZT-PEMFCs.  In  this  study,  a  novel  design 
of  a  micro-diaphragm  flow  channel  with  piezoelectric  device  for 
PEMFCs  is  proposed.  A  transitional  three-dimensional  model 
[31]  is  presented  to  calculate  the  performance  of  PEMFCs  at 
different  PZT  frequencies. 

2.  Mechanisms  of  PZT-PEMFC 

The  actuating  micro-diaphragm  with  piezoelectric  effects  is 
utilized  as  one  of  the  flow-channel  PEMFC  systems,  termed 
the  PZT-PEMFC,  shown  in  Fig.  1.  This  newly  designed  gas 
pump  with  a  piezoelectric  actuation  structure  can  feed  air  into 
the  system  of  an  air-breathing  PEMFC.  While  the  actuator  is 
moving  outward  to  increase  the  cathode  channel  volume,  the  air 
is  sucked  into  the  chamber;  when  moving  inward  to  decrease 
the  cathode  channel  volume,  the  air  is  compressed  into  the  cat¬ 
alyst  layer  and  the  chemical  reaction  is  enhanced,  as  shown  in 
Fig.  2. 


404 


H.K.  Ma  et  at.  /  Journal  of  Power  Sources  180  (2008)  402M-09 


MEA 


Valve:  Open 


—I 


ft 

— .x//^ 

1 

m 

ft  t  t  t  t 

■J _ 

%  II 

Valve:  Close 


(b)  Rejected  gases  out  of  the  channel 
Fig.  2.  Actuating  mechanisms  in  a  PZT-PEMFC  system. 


2.1.  Air-breathing  PZT-PEMFC  cycle 

The  air-standard  PZT-PEMFC  cycle,  shown  in  Fig.  3,  is  an 
ideal  cycle  that  describes  the  processes  of  an  air-breathing  PZT- 
PEMFC.  Process  1-2  is  an  expansion  process  as  the  PZT  moves 
outward.  The  pressure  in  the  cathode  channel  is  lower  than  the 
atmospheric  pressure,  P$.  In  the  constant  volume  process  2-3, 
air  is  sucked  into  the  cathode  channel  at  a  lower  pressure.  Process 
3-4  is  a  compression  process  where  the  PZT  moves  inward  and 
the  increasing  pressure  compresses  more  oxygen  into  the  cata¬ 
lyst  layer.  In  the  constant  volume  process  4-1,  water  vapor  and 
un-reacted  reactants  are  pumped  out  from  the  cathode  channel. 


2.2.  Actuation  analysis  of  PZT-PEMFC 


The  driving  forces  at  the  anode  and  cathode  inlets  are  dif¬ 
ferent.  Inlet  hydrogen  velocity  at  the  anode  is  assumed  to  be 
constant.  On  the  other  hand,  inlet  air  velocity  at  the  cathode 
is  driven  by  PZT  vibrations.  The  equation  of  PZT  motion  is 
assumed  to  be  the  step  function  shown  in  the  following  equation: 


VpZX  =  —  (-0.0005  x  STEP 
d  t  l 


sin 


( 27 root 


?)]} 


(i) 


o 

— 


a. 

E 

< 


Inflow 


Outflow 

Time 

Fig.  4.  The  step  function  of  periodical  flow  at  cathode. 


Also,  the  inflow  and  outflow  periods  in  the  channel  induced  by 
the  step  function  are  shown  in  Fig.  4.  The  inlet  air-flow  rate  at 
the  cathode  is  expressed  by  the  Reynolds  transport  theorem  as 


DM 

~Dt 


Dt 


=  l[  PdV+  [ 

ot  Jew  Jc 


pVpzj'n  d A 


CS,PZT 


+ 


[  pcVcmndA  =  0 
J  CS,in 


(2) 


The  gas  density  is  expressed  by  the  equation  of  state  for  an  ideal 
gas: 


(3) 


Substituting  Eq.  (3)  into  Eq.  (2),  the  mass  conservation  equation 
is  derived  as 


1  8  f 

R  dt  Jew 


p 

T 


dV  +  Pcs,pztEpztApzt 


PcVcAm  —  0 


(4) 


The  inlet  velocity  at  the  cathode,  Vc,  is  derived  in  the  following 
equation: 


'1  dr 

R  dt  Jew 


P 

— dV  +  pcs  kpzT  Apzt 


(5) 


3.  Theoretical  model  for  PZT-PEMFC 


A  numerical  model  of  the  novel  fuel  cell,  based  upon  the 
SIMPLEC  procedure,  has  also  been  developed.  The  assumptions 
are  as  follows: 


(1)  The  reactants  and  products  are  treated  as  ideal  gases. 

(2)  The  Stefan-Max  well  equations  are  applied  to  multi- species 
diffusion. 

(3)  Contacting  electric  and  heat  resistances  among  the  chan¬ 
nel  layer,  diffusion  layer,  catalyst  layer,  and  membrane  are 
ignored. 

(4)  The  Nernst-Planck  equation  is  used  for  the  proton  transport 
through  the  membrane. 

(5)  Ohm’s  law  is  applied  across  the  entire  region  of  the  polar¬ 
ization  curve. 

(6)  The  porosity  and  permeability  of  the  porous  media  are  uni¬ 
form,  and  the  porous  media  is  isotropic  and  homogeneous. 

(7)  The  effect  of  gravity  is  ignored. 


Fig.  3.  Pressure- volume  diagram  of  PZT-PEMFC  cycle. 
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PZT:  0.1mm 

Cathode  Collector:  1.0  mm 


MEA:  0.573  mm 


Anode  Collector:  1.5  mm 


Fig.  5.  Dimensions  of  the  air-breathing  PZT-PEMFC. 


Table  1 

Gas  composition  of  fuel  and  air  at  inlet  conditions 


Fuel 

Hydrogen 

73% 

Water  vapor 

27% 

Nitrogen 

70% 

Air 

Oxygen 

25% 

Water  vapor 

5% 

Unit:  molar  concentration. 


(8)  The  amplitude  of  the  PZT  device  is  assumed  as  a  constant 
at  different  frequencies. 

(9)  The  water  existing  in  the  fuel  cell  is  assumed  to  be  100% 
water  vapor. 

The  parameters  of  MEA  refer  to  a  DuPont  membrane,  and  the 
dimensions  of  the  PEMFC  are  shown  in  Fig.  5.  The  composition 
of  the  inlet  humidified  fuel  and  air  are  shown  in  Table  1 .  In  addi¬ 
tion,  the  inlet  fuel  velocity  is  in  the  range  of  0.25-1  ms-1.  The 
inlet  air  velocities  vary  with  different  PZT  frequencies,  which 
are  shown  in  Table  2. 

3.1.  Channel  layer 

The  governing  equations  in  the  channel  include  the  continu¬ 
ity,  momentum,  and  species  equations.  Since  the  driving  forces 
of  the  inlet  gas  are  different  at  the  anode  and  cathode,  the  conti¬ 
nuity,  momentum  and  energy  equations  should  also  be  separated 
into  two  sides: 


where  p  is  the  mixture  fluid  density,  which  can  be  determined 
by  the  following  equation: 

N 

P  =  5>  (7) 

1=1 

The  momentum  equation  is  written  as 

UpV a)  +  (VpVaVa)  =  ~Vp  +  Vr  (8) 

ot 

The  energy  equation  is: 

|(/oA)  +  V*(pVaA)  =  V*?  +  t:  VVa+^  (9) 

3 1  d  t 

where  q  can  be  expressed  by  the  following  equation: 
q  =  keftVT  +  J2  J,h,  (10) 

i= gas 

where  i!  Her  is  the  generated  heat  which  is  induced  by  the  current 
flow. 

The  species  equation  is: 

^(pF,)  +  V(pVa^)  =  VVf  (11) 

ot 

where  F;  is  the  mass  fractions  in  the  gas  phase  of  the  ith  species. 
Ji  is  the  diffusive  flux  which  can  be  expressed  by  the  following 
equation: 

Ji  =  pA'.effVy,  +  DLeffV M  -  pM^DuaVYj 

j 

-  pVM^D^ffVF;  (12) 

j 

where  pA;effVF*  is  the  Fickian  diffusion  coefficient,  which  is 
induced  by  concentration  difference.  Other  terms  are  the  cor¬ 
rections  for  multiple  species  within  the  porous  medium.  The 
fluid  diffusion  coefficient  should  be  modified  to  which  is 
derived  by  Bruggeman.  [32] 


3.1.1.  Anode 

The  inlet  velocity  of  the  anode  side  is  constant,  therefore,  the 
continuity  equation  can  be  written  as 

J  +  V(pVa)  =  0  (6) 

ot 


Table  2 

Operating  parameters 


Anode 

inlet 


Cathode 

inlet 


Fa 

Pressure 

Temperature 

Fc 

Pressure 

Temperature 


0.25  ms *  1 
1  atm 
353  K 

Eq.  (8),  ms-1 

1  atm 
353  K 


3.1.2.  Cathode 

Most  governing  equations  for  the  cathode  channel  are  similar 
to  the  ones  for  the  anode  channel.  However,  the  driving  force  of 
the  cathode  gas  is  given  by  the  PZT.  Therefore,  the  continuity, 
momentum,  energy  and  species  equations  in  the  cathode  channel 
should  be  modified  as  follows: 


3p 

f  +  VOoVc)  =  o 

ot 

UpV, c)  +  (VVVcVc)  =  -Vp  +  v*r 

ot 

U ph )  +  V(pVcA)  =  V*  +  r  :  Wc  +  ^ 
3 1  at 

l(pYi)  +  V'(pVcYi)  =  Wi 
at 


(13) 

(14) 
05) 
(16) 
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3.2.  Diffusion  layer 


The  species  equation  is: 


For  porous  media,  the  governing  equations  should  be  mod¬ 
ified  by  the  porosity,  s,  and  permeability,  k  of  the  media.  The 
governing  equations  are  as  follows: 

3.2.1.  Anode 

3. 2. 1.1.  Continuity  equation.  The  inlet  velocity  of  the  anode 
side  is  constant,  therefore,  the  continuity  equation  can  be  written 
as 


-(ep)  +  V(epVk)  =  0 

ot 

The  momentum  equation  is: 


(17) 


l(spV.A)  +  (VVpVaVa)  =  -eVp  +  V*(er)  +  (18) 

at  k 

where  s2  p,V  k-1  is  Darcy’s  drag  force  describing  the  fluid  flow¬ 
ing  through  the  porous  media.  The  mixture  viscosity  is  /x,  which 
can  be  determined  by  the  following  equation 


N 

r-  = 

i= 1 


(19) 


The  species  equation  is: 


-(epYi)  +  V(epVfiYi)  =  VJi 
ot 


(21) 


3.2.2.  Cathode 

The  cathode  inlet  velocity  is  affected  by  the  PZT  device. 
Thus,  the  continuity,  momentum,  energy,  and  species  equations 
are  shown  as 


—  (ep)  +  V(epVc)  =  0 

Ot 

2(spVc)  +  (VspVcVc)  =  -eVp  +  V*(er) 

Ot 

l~(eph)+V(epVch)=Vq+er  :  VVc+e^ 
3 1  at 


l-ispYi)  +  ViepVYi)  =  VV,  +  oj, 

Ot 


(26) 


where  coi  is  the  production  rate  during  the  gas  phase  of  the  ith 
species. 

In  electrochemical  reactions,  the  production  rate  of  each 
species  in  the  gas  phase  can  be  expressed  as 


COi 


i  -  «  -  Of 


(27) 


where  at  is  the  concentration  exponent  and  the  transfer  current 
/T  can  be  expressed  by  Butler- Volmer  equation  as  in  Eq.  (28) 

n  ak 

k=  1 

(28) 


JT  = 


JO 


nt|l^,ef] 


oik 


exp 


OigF 

RT 


r\  I  —exp 


acF 
- 1 

RT 


where  r)  is  an  over  potential,  which  is  the  difference  between  the 
solid  potential  cps  and  the  porous  (or  liquid)  potential  <^f- 
A  set  of  non-linear  equations  for  mass-fractions  at  the  pore- 
catalyst  interface  can  be  obtained  by  substituting  Eqs.  (27)  and 
(28)  into  the  following  equation: 


where  sub /  is  the  fluid  species,  and  7/  and  /x*  are  the  mass  fraction 
and  viscosity  of  the  ith  species,  respectively. 

The  energy  equation  is  written  as 

3  -+  dp  imi 

—(sph)  +  V(epVah)  =  Vq  +  ex  :  V Va  +  e- - jTq  +  — 

dt  at  a 

(20) 


u>i  =  pDj 


YfJ  ~  Yc 


VJe  ff 


(29) 


Since  the  volume  of  the  catalyst  layer  is  much  smaller  than  the 
entire  fuel  cell,  the  effective  surface-to-volume  ratio,  [A/V]e ff, 
should  be  considered.  However,  the  gas -production  rate,  &>;,  is 
only  applicable  in  the  catalyst  layer  and  when  the  other  oofs  are 
0. 

The  current  conservation  equation  is  given  by: 


Vi  =  0 


(30) 


Since  the  current  flows  through  the  porous  medium,  the  current 
can  be  divided  into  two  parts: 


i  =  if  +  is 


(31) 


where  if  is  the  current  flowing  through  the  fluid  and  is  is  the  cur¬ 
rent  flowing  through  the  solid.  Since  they  are  porous  mediums, 
the  electrons  transfer  from  fluid  to  solid  and  solid  to  fluid.  Thus, 


-(epYi)  +  V(epVcYi)  =  V/f 

ot 

3.3.  Catalyst  layer 


The  important  governing  equations  in  the  catalyst  layer  are 
the  species  and  current  conservation  equations.  Since  the  chem¬ 
ical  reaction  and  current  generation  take  place  in  this  layer,  the 
species  and  current  conservation  equations  should  be  written  as 
below. 


(22) 

the  transfer  current  can  be  written  in  the  following  equation: 

e2ptVc 

k 

(23) 

h  =  Vi*  =  -Vif 

From  Ohm’s  law,  Eq.  (32)  can  also  be  changed  to: 

(32) 

i*i 

~jT0+  — 

(7 

(24) 

V(crfV^f)  =  —  V*(ersV^>s)  =  jT 

(33) 

(25) 

4.  Results  and  discussion 

A  novel  design  of  PZT-PEMFCs  has  been  developed  and 
a  three-dimensional,  transitional  model  has  been  successfully 
built  to  account  for  its  major  phenomena  and  performance. 

4.1.  Pressure  field  analysis 

The  pressure  distribution  of  the  cathode  channel  in  an  air- 
breathing  PZT-PEMFC  is  shown  in  Fig.  6.  The  expansion 
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P  -  N/m*2 
90 


Fig.  6.  Pressure  distribution  of  the  cathode  channel  in  the  PZT-PEMFC. 


pressure  in  process  1-2  reaches  (Pq  —  80)  N  m-2  at  the  cathode 
channel.  In  process  2-3,  air  is  sucked  into  the  cathode  channel 
until  the  pressure  increases  to  Po .  In  the  compression  process 
3-4,  the  pressure  reaches  (Po  +  90)  Nm“2,  which  compresses 
more  oxygen  into  the  catalyst  layer,  and  thus  enhances  the  chem¬ 
ical  reaction  rate.  In  process  4-1,  the  produced  water  vapor  and 
un-reacted  reactants  are  pumped  out  until  the  channel  pressure 
is  reduced  to  Po. 

4.2.  Piezoelectric  effect  on  oxygen  profile 

At  high  PZT  frequencies  (/=  8,  16  Hz),  oxygen  input  is  much 
higher  than  at  low  PZT  frequencies  (f=  1,  2  Hz).  In  general, 
the  higher  the  frequency,  the  more  air  is  sucked  into  the  cath¬ 
ode  channel.  High-consumption  rates  (98%  and  87%)  of  oxygen 
occur  at  low  frequencies  (f=  1,  2  Hz)  owing  to  a  deficient  oxy¬ 
gen  input,  as  shown  in  Fig.  7.  In  contrast,  low-consumption  rates 
(47%  and  29%)  take  place  at  high  frequencies  (f=  8,  16  Hz)  due 
to  an  excess  supply  of  oxygen.  Thus,  a  higher  PZT  frequency 
( f-  16  Hz)  causes  a  higher  concentration  of  oxygen  at  the  cath¬ 
ode  channel,  as  is  shown  in  Fig.  8.  This  means  that  the  oxygen 
at  the  cathode  is  rich  and  actually  leads  to  higher  hydrogen 
consumption  and  better  fuel  cell  performance. 

4.3.  Piezoelectric  effect  on  water  vapor  profile 

At  low  frequencies  (f=  1,  2  Hz),  PZT  actuation  leads  to  cur¬ 
rent  output  vibration,  drained  water  vapor  fluctuation,  and  low 
hydrogen  consumption.  At  f-  1  Hz,  most  of  the  water  vapor 
drains  out  during  the  first  compression  due  to  the  initial  condi- 


Fig.  7.  Oxygen  consumption  profiles  under  different  PZT  frequencies  (0.7  V). 

tions,  as  shown  in  Fig.  9(a).  However,  Fig.  9(b)  shows  a  different 
trend  at  f-  2  Hz,  where  the  maximum  PZT  compression  and  the 
maximum  current  generation  occur  at  the  same  time.  Mean¬ 
while,  the  water  vapor  produced  in  the  first  period  cannot  be 
completely  drained,  so  the  remainder  of  the  produced  water 
vapor  is  left  for  the  second  period,  which  appears  to  have  a 
higher  gas  output. 

At  high  PZT  frequencies  (f-  8,  16  Hz),  more  water  vapor  is 
produced  due  to  an  excess  in  oxygen  input.  However,  a  major 
portion  of  the  produced  water  vapor  can  be  drained  out  of  the 
cathode  channel  by  PZT  actuation,  which  immediately  leads  to 
a  low  water  vapor  concentration  at  the  cathode,  as  shown  in 
Fig.  10. 


Process  1  -  2 


02 

0.25 


Process  2-3  Process  3-4  Process  4  -  1 

Fig.  8.  Oxygen  field  at /=  16  Hz. 
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Period  (sec.  /cycle) 

(a)/=  1  Hz 


Period  (sec.  /cycle) 

(b)/=2Hz 

Fig.  9.  The  current  and  mass  flow  rate  of  water  vapor  at  cathode  inlet  and  outlet 
(0.7  V). 


Fig.  12.  Current  profiles  with  period  at  different  PZT  frequencies  (0.7  V). 

4.4.  The  current  density  at  different  PZT  frequencies 

From  a  polarization  curve  diagram  at  low  frequency,  as 
shown  in  Fig.  11,  most  current  losses  occur  in  the  region  of 
concentration  losses  due  to  insufficient  oxygen  at  the  cathode 
channel.  Therefore,  a  higher  frequency  of  PZT  has  to  supply 
sufficient  oxygen  to  the  cathode  channel.  In  addition,  high  PZT 
frequencies  generate  a  stable  current,  which  overcome  the  con¬ 
centration  losses  and  eliminate  current  fluctuation.  Hence,  the 
PZT-PEMFC  generates  higher  current  density  at /=  16  Hz. 

For  high  PZT  frequencies,  the  limited  current  is  0.35  A  at 
Va  =  0.25ms-1,  and  the  hydrogen  consumption  rate  is  high, 
which  leads  to  concentration  losses  and  hydrogen  deficiency; 


Fig.  10.  The  water  vapor  field  at /=  16  Hz. 


(a)  f  =  2 


(b)  f  =  8 


(c)f=  16 


Fig.  11.  Polarization  (/-V)  curves  under  different  PZT  frequencies. 
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that  is,  hydrogen  input  increment  can  alleviate  the  concentration 
losses. 

Although  the  PZT  frequencies  at /=  16,  32  Hz  are  high 
enough  to  supply  more  oxygen  to  the  cathode  channel,  the  cur¬ 
rent  generation  is  limited  by  the  hydrogen  supply.  In  Fig.  12, 
the  optimal  operation  condition,  in  terms  of  maximum  current, 
is  0.17  A  at /=  16  Hz,  under  0.7  V,  with  a  hydrogen  inlet  veloc¬ 
ity  of  0.25  ms-1.  Since  the  hydrogen  is  insufficient  at /=  32  Hz, 
increasing  the  PZT  frequency  fails  to  enhance  the  performance 
any  further. 

5.  Conclusions 

A  novel  design  of  PZT-PEMFCs  has  been  developed  and 
a  three-dimensional,  transitional  model  has  been  successfully 
built  to  account  for  its  major  phenomena  and  fuel  cell  perfor¬ 
mance.  The  major  findings  are  as  follows: 

1.  The  air-standard  PZT-PEMFC  cycle  is  proposed  to  describe 
an  air-breathing  PZT-PEMFC. 

2.  The  PZT-PEMFC,  operating  in  hydrogen-lean  conditions, 
compresses  more  oxygen  into  the  catalyst  layer,  and  thus 
enhances  the  chemical  reaction  resulting  in  a  higher  current 
output.  At  the  same  time,  produced  water  vapor  is  pumped  out 
from  the  cathode  channel  during  the  compression  process. 

3.  At  low  frequencies  (f<  8  Hz),  PZT  actuation  leads  to  current 
output  vibration,  drained  water  vapor  fluctuation,  and  low 
hydrogen  consumption. 

4.  At  high  frequencies  (f>  8  Hz),  PZT  actuation  leads  to  a  more 
stable  current  output,  drained  water,  sucked  oxygen,  higher 
hydrogen  consumption,  and  also  overcomes  concentration 
losses. 

5.  The  optimal  operation  condition,  in  terms  of  the  maximum 
current,  is  0. 17  A  at /=  16  Hz,  under  0.7  V,  and  with  a  hydro¬ 
gen  inlet  velocity  of  0.25  ms-1. 
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